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ABSTRACT 

We  study  the  scattering  of  an  incident  wave  Ly  a  flexible  panel.  The  panel  vibration 
is  governed  by  the  nonlinear  plate  equations  while  the  loading  on  th('  panel,  which  is  the 
pressure  difference  across  the  panel,  depends  on  the  reflected  and  transmitted  waves.  Two 
models  are  used  to  calculate  this  structural-acoustic  interaction  problem.  One  solves  the 
three  dimensional  nonlinear  Euler  equations  for  the  flow-flehl  (uupled  with  the  plctte  eciua- 
tions  (the  fully  coupled  model).  The  second  uses  the  linear  wave  equation  for  the'  acoustic 
field  and  exprt'sses  the  load  as  a  double  integral  involving  the  panel  oscillation  (the  decoupled 
model).  The  panel  oscillation  governed  by  a  system  of  integro-differrMitial  equations  is  solverl 
numerically  and  the  acoustic  field  is  then  defined  by  an  ex))iicit  formula.  Numerical  results 
are  obtainerl  using  the  two  models  for  linear  and  nonlinear  paind  vibrations.  The  pnMlictir)ns 
given  by  the.,e  two  models  are  in  good  agreement  but  the  computational  tinu'  needed  for  tin' 
“fully  coupled  model"  is  60  times  longer  than  that  for  "the  decouph'd  moilel". 
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1.  Introduction 


We  consider  the  scattering  of  an  incident  acoustic  wave  by  a  planar  interface,  the  (.r. 
z)  plane.  The  interface  is  a  rigid  surface  except  on  the  flexible  thin  panel.  V,  with  i>anel 
thickness  h  much  smaller  than  the  size  L  of  the  panel,  or  the  size  of  T>, 

h/LCl.  (1.1) 

This  is  a  basic  structural-acoustic  interaction  problem  simulating  the  transmission  of  an 
external  acoustic  source,  say  the  engine  noise,  through  the  airframe  into  its  interior.  It 
is  an  interaction  problem  because  the  incident  wave  excites  panel  oscillation  which  in 
turn  induces  transmitted  and  diffracted  waves.  Meanwhile,  these  waves  contrilnite  to  the 
loading  on  the  panel.  Although  in  an  engineering  problem,  the  typical  amplitude  of  the 
transverse  displacement,  r/*,  has  to  be  much  smaller  than  the  size  of  the  panel,  i.  e., 

r/7L<l,.  (1.2) 

the  problem  can  be  nonlinear  in  three  different  aspects: 

(1)  The  panel  oscillation  can  be  nonlinear  when  the  transverse  displacement  of  the  panel 
is  not  much  smaller  than  the  panel  thickness  /?,  that  is 

(1.3) 


Nonlinear  plate  equations  are  needed  usually  when  the  incident  wave  is  in  resonance 
with  the  panel  oscillation. 


(2) 


(3) 


In  the  near-field,  of  length  scale  L,  the  flow-field  is  linear  and  ol)eys  tlu'  .sinii)le  wave 
eqviation  but  in  the  far-field,  of  length  scale  nmch  larger  than  I.  the  flow-field  may 
become  nonlinear  when  the  second  order  terms  are  needed  to  account  for  the  gradual 
steepening  of  compression  waves. 

In  case  the  flow-field  is  nonuniform  and/or  the  initial  pressure  variation  is  no  longer 
much  smaller  than  the  ambient  pressure  then  we  have  no;di?iear  flow  in  tli<  neai -field. 
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We  have  to  solve  the  Euler  equations  with  appropriate  initial  data  t(»  simulate  an 

incident  wave  and  far-held  conditions. 

Regardless  of  whether  the  flow-held  is  linear  or  nonlinear,  the  plat('  should  have  sufficient 
rigidity  so  that  the  transverse  displacement  of  the  plate  remains  small  and  hence  the 
governing  ecjuations  for  the  plate  derived  under  assumption  (1.1)  remain  valifl. 

For  the  problem  of  sound  transmission  through  a  moving  airframe,  the  external  acous¬ 
tic  field  is  governed  by  the  convective  wave  equation  while  that  in  the  interior  liy  the  simple 
wave  equation.  Also  the  pressure  and  density  in  the  exterior  can  be  much  lower  than  those 
in  the  interior.  In  the  simple  model  simulating  the  experimental  data  with  no-flow  on  both 
sides  of  the  panel,  both  acoustic  fields  are  governed  by  the  simple  wave  equation  with  the 
same  ambient  pressure  and  density.  The  validity  of  the  linear  theory  for  the  flow-field  is 
confirmed  by  the  investigations  of  French  et  al.  [1,2]  and  Maestrello  et  al.  [3].  They  showed 
that  even  when  the  panel  vibration  is  nonlinear,  the  acoustic  field  can  remain  linear. 

The  pressure  difference  across  the  j^anel  induces  the  panel  oscillation  which  in  turn 
excites  the  transmitted  wave  and  an  additional  reflected  wave.  Thus  the  solution  of  the 
panel  oscillation  is  coupled  with  the  solutions  of  the  scattered  and  transmitted  waves  which 
satisfy  the  far-fleld  radiation  conditions  for  outgoing  plane  waves  (see  for  example  [4]).  In 
general  numerical  solution  of  the  scattering  of  an  incident  wave  by  an  elastic  scatterer, 
an  interface,  or  a  panel,  requires  the  introduction  of  a  finite  computational  domain  for 
both  acoustic  fields.  Higher  order  radiation  conditions  were  derived  [5-6]  so  that  they  can 
l>e  imposed  on  the  boundary  of  the  finite  computational  domain  to  give  more  accurate 
approximation  to  the  solution  in  the  unbounded  domain.  Since  the  size  of  the  computa¬ 
tional  domain  has  to  Ije  much  larger  than  the  size  of  the  panel  or  a  scatterer,  numerical 
solution  of  this  three-dimensional  unsteady  problem  is  very  teflious  especially  wIkui  the 
computation  lias  lo  lic  continued  for  a  long  time  relative  to  tin'  pi'iiod  of  oscillation  of  tlu’ 
panel.  Furthermore,  the  accuracy  of  the  solution  depends  not  fudy  on  that  of  the  iiunKui- 
cal  solution  of  the  differential  equations  but  also  on  the  aiiproximati*  Ixmndary  conditions. 
Refinements  of  the  grid  size  and  time  step  li.'.ve  to  witli  an  enlaror,Mnent  of  the  compma- 
tional  domain.  It  is  desirable  to  find  exact  boundary  conditions  so  that  the  conquitational 
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domain  does  not  have  to  be  mneli  larger  than  the  seatterer  and  is  independent  of  tlie  rlioire 
of  grid  size.  The  exact  conditions  were  presented  in  [7]  for  a  seatterer  or  inhomogeneous 
medium  of  bounded  support.  With  the  seatterer  inside  the  computational  domain,  tlie 
integral  representation  of  the  solution  of  the  simple  wave  equation  is  applieal)le  to  the 
region  outside  and  on  the  lioundary  of  the  computational  domain. 

Here  the  seatterer  is  a  panel  embedded  in  a  rigid  plane.  The  integral  representation 
can  be  applied  on  the  panel  and  the  integrand  involves  only  the  normal  velocity  of  the 
panel  [8,  9].  Thus  we  have  a  system  of  integro-differential  equations  for  the  panel  oscillation 
de<'oupled  from  the  acoustic  field.  This  .system  is  referred  to  as  the  decoupled  model. 

For  the  case  that  the  acoustic  wave  length  is  much  smaller  thaji  X,  the  integral  repre¬ 
sentation  can  be  approximated  by  the  local  derivatives  of  the  panel  oscillation  and  hence 
the  integro-differential  ecpiations  are  reduced  to  jjartial  differential  equations  [8,  9].  This 
approximation  is  not  applicable  here  because  iir  the  experimental  data,  the  acoustic  wav'e 
length  is  comparable  to  L.  We  need  numerical  solution  of  the  integro-differential  equations 
for  the  panel  oscillation  and  then  we  can  use  the  integral  representation  to  evaluate  the 
acoustic  field  and  the  nonlinear  waves  in  the  far  field  by  Whitham’s  theory  or  matched 
asymptotics. 

Note  that  the  solution  of  the  decoupled  panel  oscillation  is  an  unsteady  two  dimen¬ 
sional  problem  in  the  finite  domain  P  and  is  .several  orders  of  magnitude  simpler  than  tlu' 
solution  of  the  fully  coupled  problem.  To  show  the  efficiency  of  the  decoupled  model  and 
the  accuracy  of  the  solution  even  when  the  panel  oscillation  is  nonlinear,  we  comi)are  the 
solution  of  this  decoupled  system  with  the  solution  of  the  fully  coupled  ])rol)lem  for  which 
we  solve  the  three-dimensional  nonlinear  Euler  equations  for  the  flow-field  cou])led  with 
the  nonlinear  equations  for  the  jianel  oscillations. 

In  the  next  section  we  inesent  the  mathematical  formulations  of  these  two  models  and 
discuss  their  ajjph.  ations  in  gem'ral.  Section  3  descril)es  the  various  numerical  technitpu's 
used  to  s'i’ive  the  pit)b'eT:r'  iircolvt'fi.  1  he  results  and  discu'^sion  are  givt  u  hi  .-.t  <  iloa  4,  and 
the  conclusions  are  in  section  5. 
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2.  Formulation  of  the  Two  Models 


The  physical  problem  being  studied  experimentally  is  that  of  the  nonlinear  oscillations 
of  a  typical  aircraft  panel  excited  by  harmonic  plane  waves  at  normal  incidence.  The  panel 
is  clamped  onto  a  large  rigid  plate.  Let  the  typical  panel  be  represented  by  a  rectangular 
domain  T>  in  the  {x,z)  plane,  with  length  L,  width  W  and  one  vertex  located  at  (.ro,Co), 
i.  e., 

V  =  {(;r,c)  I  j”o  <  <  xo  +  I  ,  zo  <  z  <  zo  +  IT}  .  (2.1) 

The  compliment  of  T>  in  the  (x,^)  plane  represents  the  rigid  plate.  The  incident  side  is 
the  half  space  y  >  0  and  the  transmitted  side  is  the  half  space  y  <  0.  Both  sides  have 
the  same  ambient  pressure  pcc  and  temperature  Too  and  hence  the  same  speed  of  sound  C 
and  density  poo-  We  denote  the  pressure  and  velocity  potential  of  the  flow-field  by  p  and 
and  use  the  superscripts  and  —  to  denote  the  quantities  on  the  sides  y  >  0  and  y  <  0. 
On  the  side  y  >  0,  the  superscripts  i  and  r  are  used  to  denote  the  quantities  associated 
with  the  incident  and  reflected  waves,  respectively.  These  two  waves  are  mirror  images 
with  resjject  to  the  (x,  x)  plane.  When  the  incident  waves  are  plane  waves  advancing  in 
the  direction  opposite  to  the  y  axis  and  hitting  the  (x,  plane  at  f  =  0,  we  have 

4)^'\t,.r,y,z)  =  f{Ct  +  y)  and  x,y,  z)  =  f{Cf  ~  y)  .  (2.2) 


with  /(^)  =0  for  ^  <  0.  Therefore  =0  when  0  >  y  <  Ct. 


Transverse  oscillation,  r/,  of  the  rectangular  flexibh’  panel,  P,  is  excited  l>y  tin-  pressure 
diffe  rence  across  the  panel  for  t  >  0  and  the  oscillation  in  turn  induces  scattered  waves  o’ 
in  the  incident  side  and  transmitted  waves  Both  waves  have  th(“  homogeneous  initial 
data  at  f  =  0, 

-  0  ,  y  >  0  , 

(2.3) 

Under  linear  theory,  we  have 


<i)''''\t,x,y,z)  =  ())^'\i,x,-y,z) 
<l)^‘\t,x,-y,z)  -  y,  c)  for  y  >  0  . 


(2.4) 
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The  wave  fronts  of  and  are  the  envelopes  of  the  sonic  sphei’es  with  radius  Ct  and 
centered  on  the  panel  in  y  >  0  and  y  <  0,  respectively.  The  pressure  in  the  inciclcnt  and 
transmitted  sides  are 

p+  -  Poo  =  ,  y  >  0 

(2.5) 

p  -poo  =  ,  y  <  0 

and  the  pressure  (lifference  across  the  panel  is 

Ap  =  p~  -  p"*"  =  2p^'^  -  [p^‘^ +p^’’^]  ,  y  =  0,  {x,z)eT>.  (2.6) 


In  this  section  we  formulate  two  mathematical  models  to  compare  with  and  compli¬ 
ment  the  experimental  studies.  We  pre.sent  the  system  of  equations  governing  the  nonlinear 
flow-field  and  panel  oscillation  (the  fully  coupled  model)  in  Section  2.1  and  formulate  the 
system  of  integro-differential  equations  for  nonlinear  panel  oscillation  inchiding  the  effec’t 
of  a  linear  acoustic  field  (the  decoupled  model)  in  Section  2.2. 

2.1  The  Fully  Coupled  Model 

In  this  model  the  structural-acoustic  interaction  is  analyzed  by  solving  the  three- 
dimensional  nonlinear  Euler  equations  together  with  the  nonlinear  plate  equations.  The 
configuration  of  the  computational  domain  is  shown  in  Fig.  1.  In  cartesian  coordinates,  x, 
y  and  z,  the  comijressible,  laonlinear  Euler  equations  can  be  written  in  conservation  form 
as 

Qt=Fr+Gy  +  H,,  (2.7) 

where  Q  is  the  vector  (p,  pu,  pv,  piv,  e),  p  is  the  density,  pit,  pv  and  pie  are  the  x,  y  and 
z  momenta  respectively,  and  e  is  the  total  energy  per  unit  volume  given  liy 

e  =  ip(it^ -f  -f  le^) -f  pr„T,  (2.8) 


with  c„  being  the  specific  heat  at  constant  volume.  In  eq.  (2.7),  the  functions  F,  G  and 


H  are: 


and 


G  = 


H  - 


/  f'  N 

put) 
pain 
\n(f  +/))/ 

(  \ 

'  pun  ‘ 

4-  }) 

pvto 
\  v((  4-  p)  ) 

(  \ 

putt' 

pvu' 

.2 


p  w  4-  p 
\  ir(  r  4-  p)  / 

In  addition  to  (ui.  (2.7),  the  efination  of  state  of  an  ideal  gas  is  ns('d: 


P  =  ,RT. 


(2.9) 


(2.10) 


where  p  is  the  pressnre,  R  the  gas  eonstant.  and  T  th<'  teiii])eratnr-“.  Since  we  assnme  that 
the  incident  wav<-  liits  the  panel  at  1  —  0.  the  How  on  the  transinitt('d  sid(‘,  y  <  0,  is  at 
rest  for  /  <  0.  The  initic'il  data  at  /  =  0  for  y  <  0  ar  ■ 

a  =r  (t  =  t(t  =  0  ,  =  Pro  iiiid  p  =  Poo  ■  (2.11) 

On  the  rigid  i)hite  we  have  zero  normal  velocity, 

e(.r,0,r)  =  ()  for  (2.12) 


The  motion  of  the  He.xihle  panel,  "P,  is  descril>ed  by  a  system  of  thn'e,  nonlinear  partial 
difh'rential  equations  given  by  [10]; 


PV'S/  4-  Pphijtt  4-  7'/f  =  ^P 

Eh  ,  1  , 

I-  (a>"  +  )(//,;  +  I'If.r) 

f  (1  -  4-  fc"  4-  t/j?/;)]  (2.13) 


0 


«rjr  +  +  (hw].,  =  -t]x(Vxx  +  “  (hlUVxz 

7t>°,  +  d\wlj.  +  d'zU^^,  =  -n:{l]zz  +  d\1]xx)  -  d'^IJxVxz 


wlu'n' 


V  7/  —  f)zxx.x  ^^jxxzz  ^Izzzz 

,  l-u  ,  1  +  ;/  „ 

d]  =  — T-.  (k  =  = 


12(1  -;/2)’ 


(2.14) 

(2.15) 


(/*’  <111(1  //;'*  lire  the  in-plaiK'  disphieemeiits,  and  7/  is  th('  trausv(^rse  displacement.  The 
physical  constfints  for  the  panel  appearing  in  (2.13)  are;  the  stiffness  (D),  the  density  (/?;»), 
th('  thickness  (/i),  the  physical  damping  (7),  the  modulus  of  elasticity  (E)  and  the  Pois.son 
nitio  (;/).  The  system  of  equations  (2.13)  is  solved  subject  to  the  homogemeous  initial 
condition 

at  t  =  0  7/”  =  v)^  =  7/  =  7/(  =  0  (2.16) 

and  th('  clamped  bouiuh'iry  conditions  on  C  given  by 

X  =  .ro  ,  ,ro  +  I  ,  77”  =  7/)"  =  7/  =  7/j.  =  0  , 

Z  =  ,  Zq+W  ,  77°  =  7e”  =  7/  =  7/.  =  0  . 


(2.17) 


In  (■(].  (2.13),  the  load  Ap  defined  by  (2.6)  contJiins  the  coupling  with  the  acoustic  field  77* 
.and  the  forcing  term  (77^'*  +7J*^0  which  re]>resents  the  load  on  the  panel  if  the  panel  were 
rigid.  We  assuiiK’  the  forcing  term  to  be  <1  luirmonic  wfive  of  the  form 


77^1  -f-  H{t)  7711  7/  =  0”^  . 


(2.18) 


where  f  and  uj  are  tlu’  amplitude  find  fre7[U('ncy  of  the  wave  and  H  denotes  th('  Hefivisich' 
unit  step  function.  The  loail  on  the  panel  P  given  by  (2.6)  bf'conu’s 


A77  =  [277^  —  f  .S777(u.'t)]  H{t)  . 


(2.19) 


Another  condition  coiqiling  the  flciw-fiidds  find  the  imnel oscillation  is  the  kiimmatic  con¬ 
ditions, 

Jl,{t,x,z)  =  v^{t,.r,0^,z)  ,  ix,z)  e  V  .  (2.20) 
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They  are  imposed  on  y  =  0*  because  of  the  small  panel  displacemt'iit ,  (11).  The  numei  ical 
scheme  for  the  solution  of  the  nonlinear  plate  equation  (2.13)  and  tin*  Euler  ('(juations  (2.7) 
for  the  transmitted  waves  with  homogeneous  initial  data  will  bt'  described  in  Sec.  3 

2.2  The  Decoupled  Model 

Now  we  analyse  th«'  nonlinear  panel  oscillation  exciti'd  by  a  w^^fk  incident  pressiire 
wave,  under  the  assumption  that  the  pressure  fluctuation  remains  much  smaller  than 
the  ambient  pressure  i. 

(/' -  Poc)/p-x  <  1  •  (2.21) 

Consecpiently,  the  small  disturbance  theory  is  api)iicab]e  to  the  flow  fit'lds  in  the'  incident 
and  transmitted  sides  and  the  pan<'l/a<’oustic  intera<'tion  problt'in  is  (h'scril)ed  l)y  (2.1) 
(2.G).  The  velocity  pot<unial  $*(/,./•.  y.  r )  is  goverm'd  by  the  simi)l(  wave  ecpiation. 

-  Oi,  -  -  di, )  =  0  .  for  ±  y  >  0  .  (2.22) 

find  the  acoustic  pn'sstire  and  velocity  are  relatefl  to  the  potemtial  by 

p±  =  ami  V*  =  V‘I>*  .  (2.23) 

In  particular,  we  have  v  =  where  r  (hmotes  tlu'  vertical  V('locity,  The  vertical  dis¬ 

placement  of  the  i)anel,  ;/(/,  ,r,  z),  is  governe(l  by  the  .system  f  partial  differential  ecpiations 
(2.13).  To  produce  the  si)ecial  fencing  teTin,  (;/'*  -f-  acting  on  tin*  i)anel  specified  by 

(2.18)  in  Sec.  2.1.  the-  incident  {jotenitial  should  be 

0<‘’(f..r.y,c)  =  - - -  cos[u.'(f  -f  ^)]H(t  +  ^)  ■  (2.24) 

Note  that  the*  incident  wave  is  a  solution  of  (2.22)  in  the  whole  space.  When  the  panel 
is  rigid,  th(>  incident  wave  is  reflected  by  the  rigid  (.r,  z)  i)lane  and  the  rt'ttecteal  wa^’*'  is, 

,  s,  y.  z)  =  .  X, —y,  z)  for  y>0.  (2.25^ 
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Oil  the  side  y  <"  0.  the  flow-field  reniaiiis  at  rest  heraiise  the-re  is  no  traiisinitted  wave. 

For  a  flexible  jiaiiel,  the  panel  oscillation  excited  liy  the  jin'ssure  difference  across  the 
panel  — .r,  0“^.  c)  induces  a  scattered  wave  on  the  incident  side  and  a  transmitted 
wave  Tims  we  write 


$  =  0(‘> +  ,  y>0 


{2:2Ga) 


and 

$  =  y<0.  {2.26b) 

Note  that  ^O)  governed  by  the  simple  wave  equation  (2.22)  in  y  >  0  and  y  <  0 

respectively.  The  kinematic  conditions  on  the  (x.  c)  plane.  (2.20)  and  (2.23).  become 

.r.O'^ .  z)  =  dcT/(t.  .r,  :)  and  c)f,o''^(f,.r.0~ .  z)  —  0tf/(f..r.  z)  .  (2.27) 

where  y  denotes  the  extension  of  the  transverse  displacement  of  the  pant'l.  V.  to  the  (.r. 
c)  jilane.  i.  e.. 

7/  =  0  .  {.r.  z)  ^  72  .  (2.2S) 

If  the  incident  wave  front  hits  the  iianel  at  t  =  0.  we  can  imi)os('  the  homogeneous  initial 
conditions  on  0b'')_  <^(0  7^ 

(vf)*’*  =  0  .  dt(f)^^^  =  0  ,  0^*  =  0  .  dt<P^'^  =  0  and  y  =  0  .  Oty  =  0  for  f  <  0  . 

(2.29) 

Since  the  ambient  fluid  alxwe  and  b«-low  the  panel  are  the  srnne.  (2.27)  imiilies  that  the 
velocity  potential  induced  by  the  panel  oscillation  has  to  be  anti-symnu'tric  in  y.  i.  e.. 

0b*(#.  ,/•.  y,  ;)  =  X,  —  ij,  z)  .  y  <  0  .  (2.30) 

Then  the  pressure  difft'ience  across  can  be  written  as 

Ap  =  2p.^[5,<;6b'(#,.r,0+..')-f  c9,0b”(/.-r-O+.r)]  .  (.r.z)e  'D  .  (2.31) 
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The  (Icrivatioii  of  the  sy^tnn  of  ot[uations.  I ’2.21)  to  (2.31).  ('(laiplotcs  the  foiinulariou  of 
’’(■  stnu'tiu'aly^acoii.sric  iiitoiai’f ion  i)i'ol)lo:n  liimar  theory  for  the  a<'oustie  fields.  Tlie 

panel  oscillation,  which  can  la-  nonlinear  is  aoverneil  hy  the  syst(‘in  of  e(piati()ns.  i2.1..1)  - 
(2.17.1. 


The  velocity  jioti'ntial  C)*"*  induced  hy  the  pan<  1  oscillation  //  is  lioverned  hy  the  wave 
efiuation  (2.22).  the  initial  .  onditioiis  (2.29)  and  the  houndary  condition  (2.27i.  The 
solution  is  yiven  hy  the  Kirchlioff  huinula. 


>{).  =  )  = 


(2.32) 


where  n  —  I"' ~  <lenotes  thi' <lisranc(>  from  a  point  (./■.  //  >  0.  v  )  to 

a  source  at  (.r'.O.  r')  and  {■}  denotes  th<'  retarded  value  of.  i.  e.. 


H 

{ijiit . A).  :' )}  =  tjiit  -  O.v')  . 


(2.33) 


Tlu'  domain  <if  dependence  of  o* ''hi  •  .r. //.  c)  i-s  the  circuhir  disc  Ti  in  tin'  i.r'.:')  ])l<uie.  i. 
e. . 

H  I  n  <  Ct  or  =  (.;•'  —  ./’)*  +  (r'  —  c)"  C'~t~  —  ,(/"  .  (2.34) 

The  donifiin  of  integration  in  (2.32)  is  the  inter.s<'ction  of  Ti  and  the  jianel.  i.  e.. 


g  =  n  n  P  ■ 


(2.35) 


X('W  ve  introduce  the  polar  coordinates.  r.O  .  c<‘iit<'red  at  (./'.v).  i. 


(■  cos 


0  and 


r  sm  I 


(2.30: 


and  ( 2.32 )  hecouies. 


(f.  r.!)  -^O.v) 


R  ,  r  dr 

iW  fiiif  -  ~-.r.  d)  ] 


( 2.37  ) 


with 


f  r  cos  d. .:  t  rsind) 


:2..3Si 
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Here  R  =  (r^  +  represents  the  slant  height  of  a  circular  cone  with  vertex  'P{x,y,z), 

a  vertical  axis  and  a  base  circle  of  radius  r  in  the  (x,  z)  plane.  On  account  of  (2,28),  we 
can  extend  the  domain  of  integration  to  H  and  rewrite  (2,37)  as  an  iterated  integral, 

4,<’\t,x,y>0,z)  =  ~^  j  '  [  r  .  (2.39) 

The  integral  in  9  represents  the  contribution  of  the  sources  on  a  base  cir(de  of  radius  r. 

In  (2.31),  we  need  to  relate  the  unknown,  0"^, 2)  on  the  panel  to  7){t,x,y)  for 

{x,z)  G  V.  This  is  obtained  by  differentiating  the  above  equation  with  respect  to  t  and 
using  (2.28).  This  results  in  replacing  on  the  left-hand-side  of  (2.39)  and  gi  in  the 
integrand  by  and  gtt^  respectively. 

As  y  — >  0"^,  i?  — »  ?•  and  the  circular  disc  7i  is  bounded  l>y  the  .sonic  circle  ;•  =  Ct. 
Equation  (2.39)  yields 

1  f2re  fCl 

,z)  =  git{t  -  dr  d9  .  (2.40) 

Note  that  in  this  form  we  remove  the  kernel  1/i?  in  (2.32),  which  becomes  singular  as  y  — »  0 
and  7'  — »  0  and  shows  that  the  area  of  the  domain  of  integration  is  bounded  above  by  2nCt. 
By  using  (2.31)  and  (2,40),  (2.13)  becomes  a  system  of  integro-differrential  equations  for 
the  panel  oscillation,  r/.  The  initial  and  boundary  conditions  for  7/  are  (2.16)  and  (2.17). 
Thus  we  complete  the  formulation  of  the  decoupled  model.  The  numerical  solution  of  this 
system  will  lie  described  in  Sec.  3, 

2.3  Application  to  Noise  Transmission  in  a  Supersonic  Flight  at  High  Altitude 

The  analysis  reported  in  Sec.  2.2  is  based  on  tw('  assumptions 

(i)  The  flow-fields  above  and  lielow  the  panel  have  the  same  ambient  condition,  i.  e., 

pto  =  pto  =  PZo  anfl  =  0. 

(ii)  The  panel  is  planar  and  not  prestre.s.sed,  i.  e,,  under  zero  tension. 
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The  analysis  can  lie  employed  to  estimate  the  noise  transmitted  into  the  fuselage  flying 
at  low  speed  and  a  relatively  low  levri,  where  the  static  pressure  remains  nearly  ('cjual  to 
one  atmosi)here.  Assumption  (i)  is  valid  when  terms  of  the  order  of  the  flight  Mach  number 
M  are  omitted.  Assumption  (ii)  is  valid  since  the  size  of  the  panel  L  is  much  smaller  than 
the  radius  of  the  fuselage,  Rf.  There  is  no  hoop  stress  when  =  p'^  under  assumption 
(i).  In  practice,  theie  is  a  pressure  difference  (p^  —  pi  )-  ^'ven  for  a  low  Mach  number 
flight. 

For  an  airplane  cruising  at  subsonic  or  supersonic  speed,  it  is  necessary  to  account  for 
the  differences  in  densities  and  in  pressures  inside  and  outside  the  fuselage  and  the  moving 
media  outside  and  to  include  the  contribution  of  the  hoop  stress  to  the  panel  oscillation. 
The  formulation  of  the  decoupled  models  for  structural  acoustic  interactions  in  a  moving 
media  at  sidjsonic  and  supersonic  speeds  is  carried  out  using  the  linear  theory  for  the  flow- 
field.  The  system  of  integro-differential  equations  for  the  panel  oscillation  is  then  reduced 
to  a  system  of  differential  equations  when  the  acoustic  wave  length  is  much  smaller  than 
the  surface  wave  length  [11]. 

In  this  subsection,  we  show  that  with  minor  modifications  the  above  analysis  becomes 
applicable  to  a  supersonic  airplane  cruising  at  high  altitude.  We  make  use  of  the  fact  that 
the  cabin  pressure,  p^,  and  temijerature  and  hence  the  density  are  maintained  near  the 
amluent  level  while  the  outside  press\ire  and  density  p^  at  the  flight  altitude  H  are 
much  lower.  The  density  ratio,  is  about  1/10  or  1/20  when  H  =  60,000  or  SO.  000 

ft.  The  "ame  is  ♦^r"'^  for  the  pressure  ratio  while  the  absolute  temperature  ratio  and  hence 
the  ratio  of  the  speed  of  sound  remains  0(1),  i.  e., 

=  ptolvZo  =  0(p)  and  C+/C-=0(1).  (2.41) 

Thus,  tlu'  density  ratio  p  serves  as  a  small  ijarameter  in  th»'  following  modification  of  the 
analysis  in  Sec.  2.2. 

On  the  transmitted  .side,  y  <  0,  the  governing  equations,  (2.13)  (2.23),  (2.26b),  (2.27) 
and  (2.29).  remain  valid  when  the  superscrii)t  "  is  added  to  p.y^  and  C.  On  the  incident 
side,  y  >  0.  equations  (2.25),  (2.26a),  aiul  (2.29)  remain  the  same  while  (2.22).  (2.23) 
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and  (2.27)  have*  tci  be  modified  to  include  the  convection  terin.s.  Consequently,  (2.30)  aiul 
(2.31)  will  be  modified.  In  particular,  we  note  the  acejustic  pressure  relationship  (2.23)  in 
the  transmitted  side. 


P~  =  ^  y<0  (2.42) 

anel  that  on  the  incident  side, 

P'^  -Pt=  =  -PP~^\di  +  A/C+c>,]$  ,  y  >  0  .  (2.43) 


Using  (2.26a)  and  (2.43),  we  have 

P^"^  =  -Pp:A^i  +  y  >  0  ,  (2.44) 

Likewise,  the  kinematic  condition  on  the  y  =  O'*'  is  modified  and  (2.27)  becomes, 

d,j(l)^'''\t,x,0'^ ,z)  =  [dt  +  MC'^di]7j{t,x,z)  and  dy4>^*\t,x,0~ .z)  =  diy](t,x,z)  .  (2.45) 

With  the  convection  term  at  most  of  the  order  of  the  unsteady  term,  or  one  order  smaller 
uiifler  the  short  acoustic  wave  approximation,  we  deduce  from  (2.45) 

dy<f>^^\t,x,0-  ,z)  =  0{dy<i>^^Ht,x,^'\z))  (2.46) 

and  then  from  (2.44) 

;/’>(Ux,0+,r)  =  0{fi  p<‘>(^.^0-,c))  .  (2.47) 


By  using  (2.41)  and  (2.44),  the  load  on  the  panel,  which  is  the  pressure  difference,  becomes 


Ap  =  (Pi  -  Pj.)  +  [-2p"'  +  (p'"  -  p'”)j/f(() 

=  Pill  -  0(1,)]  +  |-2p'"  +  p"'(l  +  0(p)|//(() 


(2.48) 


The  constant  load  (p^  —  p^)  »  p^  is  present  for  all  t.  The  unsteady  load  due  to  the 
acoustic  waves  is  present  only  for  f  >  0  and  is  much  smaller  then  p^.  The  constant  load 
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produces  a  static  dt'flectioii  ;)  which  does  not  generate'  far-iield  sound.  The  uiisre-ady 

load  t'xcites  panel  oscillation  7j'{t,j\z)  and  geu»‘rat<‘s  the  far-tield  sound.  Then-fore.  we 
separate  the  static  panel  dt'flt'ction  from  the  unsteady  oscillation, 

7/  =  ?/o(.c.r)  + -  (2.49) 

Note  that  the  panel  oscillation.  ?/,  is  coupled  with  only  the  acoustic  fit*ld  (2.48)  in  the 
transmitted  side,  since  p*'*  is  prescribc'd  and  p*’'  is  0(fi)  related  to  p*''.  This  is  equivalent 
to  say  that  in  a  high  altitude  flight  tlu'  acotistic  damping  in  the  flow-fie'hl  outside  of  the 
fuselage  is  of  the  ortler  of  the  density  ratio  ft  =  P^/pT  I't'lative  to  that  insirle  and  hence 
is  neglegihle. 

Since  the  fluid  is  at  rest  inside  the  ftiselage,  the  transmitted  wave  obeys  tin*  simph' 
wave  equation  (2.22).  We  can  make  use  of  the  analysis  in  Sec.  2.2,  in  particular.  (2.30) 
and  (2.40),  to  relate  the  on  surface  transmitted  pressure  to  the  panel  oscillation. 

p(')(/,.r.0-,c)  =  -p-<9,d^''  =  -^^  '  g„(t--^.r.0)dr,ie.  (2.50) 

Thus  we  have  a  closed  system  for  tin*  panel  oscillation,  r/h  from  which  we  can  then  comj)ute 
the  transmitted  and  scatt<'red  acoustic  fi<'lds. 

3.  Numerical  Methods 

For  the  ftdly  coupltd  model,  the  nonlinear  Euler  ecpiations,  (2.7).  are  solved  using  an 
explicit  finite  difference  scheiiK'.  The  .s<'hem('.  which  is  a  generalization  of  MacCormack's 
scheme  obtained  by  Gottlieb  and  Turkel  [12],  is  fourth-order  accurate  in  si)ace  and  second- 
ord<T  accurate'  in  time.  Further  eletails  on  tin*  implementation  of  the  scheme'  e-an  be'  fe)unel 
in  Freneli  ct  al.  [1,2].  The'  physie-al  boundary  of  the  e-omi)ntational  de)main  (the'  be)tte)m 
boundary,  see  Fig.  1)  is  e'ompe)seel  e)f  a  flexible  panel  e'lampe'd  be'twe'e'ii  rigiel  jelates.  Ove'i 
the  rigid  i)lates,  the'  vt'rtie-al  v«'le)city  is  zero  (e  =  0)  anel  the  surfae-e'  te-mpe-rature-  is  T,,.. 
whie'h  in  this  pape'r  is  the  same  as  the  fluiel  te'inperature'.  T^.  The  .r  anel  r  cejuqieene'nts 
e)f  the'  velocity  (ii  anel  i/i.  re'spe'ctive'ly )  are'  obtaine'el  through  line-ar  e'Xtrape)latie)n  fre)m  the- 
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interior  of  the  computational  domain.  Tlie  pressure  Ixnindary  condition  is  ol)tained  usint;; 
the  normal  momentum  equation  !>y  simply  imposing  the  normal  gradient  of  the  sum  of 
l)ressure  and  vertical  momentum  flux  to  he  zero.  i.  e.,  dy(j)  +  pv^]  =  0.  ()v(‘r  the  flexihle 
panel,  the  vertical  velocitj  is  set  to  he  eqtial  to  that  of  the  panel  and  the  temperature  is 
T,t,.  The  X  and  r  velocity  components  and  the  pressure  are  extrapolated  from  the  interior 
of  the  computational  domain.  The  ap])ropriat('  houndary  conditions  on  the  computational 
domain  are  derived  using  the  method  of  characteristics  [13].  One  should  mention  that 
the  characteristic  l)oundary  conditions  and  the  extrapolation  are  only  first  order  accurate, 
while  the  interior  scheme  is  fourth  order.  These  houndaries  are  htdieved  to  he  a  source  for 
numerical  error  in  the  fully  coupled  model.  The  dimensions  of  the  computational  domain 
are,  G1  cm,  3.05  m  and  40.64  cm  in  the  x.  y  and  r  directions,  n'sjjectively.  The  numher  of 
computational  points  used  are,  121,  241,  and  81  in  the  respective  directions. 

The  nonlinear  plate  equations,  (2.13),  ar«*  solved  using  a  finite  element  method  devel¬ 
oped  hy  Rohinson  [14].  The  panel  is  30.5  cm  long,  20.32  cm  wide  and  0.102  cm  thick,  and 
the  numher  of  elements  used  are  G  and  8  respectively. 

Since  the  grid  used  for  solving  the  jdate  equations  is  rectangular,  it  is  easlier  to 
evaluate  the  integral  in  (2.40)  in  cartesian  coordinates.  In  this  case  there  is  a  singularity 
at  each  computational  point  corresponding  to  r  =  0.  In  order  to  overcome  this  difficulty, 
a  Taylor  series  expansion  up  to  the  .second  order  is  used  to  calculate  the  contrihution  of 
the  singular  point.  The  contrihution  of  tlie  various  points  on  the  panel  that  lie  within 
the  sonic  circle  {R  <  Ct)  is  calcidated  hy  first  integrating  in  x  using  a  comhination  of 
Simpson  and  trapezoidal  rules.  The  result  is  then  integrated  in  c  using  Simi)son's  rule 
of  integration.  Because  of  the  pres«mce  of  the  retarded  tiim^  in  the  integral,  the  vertical 
velocity  of  the  plate  (?/()  is  stored  at  each  point  for  .several  time-steps.  The  evaluation 
of  the  double  integral  can  use  a  much  larger  time  step  than  the  At  for  the  integration  of 
the  differential  equations,  yet  having  the  same  degree  of  accuracy.  For  a  plate  of  givtui 
dinn'iisions  (T,  W),  and  for  a  fixed  time-step  At,  the  maximum  numh<>r  of  time-steps  to 
h('  stored  is  N  =  \/L^  +  I4^’^/(CAf ).  The  numher  N  has  to  he  changi'd  when  calculating 
the  radiated  pressure  away  from  the  plate.  The  numher  .V  is  relat('d  to  the  radius  of  the 
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largest  sonic  circle  in  the  domain  that  contains  the  flexible  panel.  For  the  cases  jjresented 
in  the  next  section  N  =  1503  for  a  time-step  At  =  4x10“*’  sec. 

4.  Results  and  Discussion 

The  numerical  schemes  for  the  two  models  presented  in  Sec.  3  are  used  to  predict  the 
vibration  of  a  flexible  panel  and  the  resulting  acoustic  radiation.  The  panel  is  forced  to 
vibrate  by  harmonic  plane  acoustic  waves  at  normal  incidence.  The  frequency  of  the  source 
is  751  Hz,  which  corresponds  to  a  natural  frequency  of  the  panel.  Two  different  amplitudes 
of  the  incident  waves  are  used.  The  properties  of  the  panel,  which  are  considered  to  be 
uniform,  are:  density  pp  =  4450.15  Kg/ni'^,  modulus  of  elasticity  E  =  1.10316x10”  N/m^, 
Poisson  ratio  =  0.33  and  a  damping  ratio  of  0.01.  The  acoustic  fluid  properties  are: 
temperature  Too  =  288.33  °K,  density  poo  =  1-23  Kg/m'^,  pressure  poo  —  1.013x10^  N/in^ 
and  sound  speed  Coo  =  340  m/sec.  The  .specific  heat  at  constant  volume  is  =  1.004 
K.J/(Kg  °K),  the  ratio  of  specific  heats  is  7  =  Cp/cv  —  1.4.  In  the-far  field,  the  fluid  is  at 
rest.  The  variables  plotted  on  figures  2-5  are  nondimensional.  The  reference  quantities  are 
given  by 


{x,y,  Z,J})ref  —  Kefi  ^ref  — 


're  f 

'c~' 


and 


Tref  = 


£! 


Pref  -  poo,  (u,V,W,^)ref  =  C,  and 


iP,^)ref  —  PooC  ■ 


(4.1) 


where  the  reference  length  is  Iref  =  0.3048  m. 


Figure  1  shows  the  configuration  of  the  computational  domain,  a  rectangular  box 
with  the  lower  side  composed  of  a  flexible  panel  clamped  between  rigid  plates.  For  a  low 
excitation  amplitude,  100  dB  or  5  x  10““*  atm..  Fig.  2  shows  that  the  panel  response  is 
linear.  In  this  case,  both  models  predict  the  same  panel  response  as  shown  by  the  figure. 
This  result  is  expected  since  the  incident  wave  is  extremely  weak.  Figure  3a-b  show  the 
time  histories  of  the  radiated  near-  and  far-field  pressures,  2.54  cm  and  1.524  in,  or  £/12 
and  5T,  away  from  the  panel  center  respectively.  The  radiated  pressure  predicted  Ijy  the 
linear  theory  (2.32)  is  in  excellent  agreement  with  that  predicted  by  the  Euler  equations 
(2.7). 
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When  the  level  of  the  excitation  is  increased  to  160  dB,  or  0.01  atm.,  the  response 
of  the  panel  becomes  nonlinear  as  shown  by  Fig.  4.  The  time  history  of  the  panel  center 
displacement  shows  a  non-periodic  behaviour  characteristic  of  nonlinear  response.  The 
figure  also  shows  that  even  when  the  panel  oscillation  is  nonlinear,  the  predictions  of  the 
two  models  are  in  reasonably  good  agreement,  because  the  pressure  variation  is  still  much 
smaller  than  poo  =  1  atm.  The  radiated  near-  and  far-held  pressure  time  histories  are 
shown  on  Fig.  5a-b.  In  the  near-held,  L/12,  the  predictions  of  the  two  models  are  in 
reasonably  good  agreement  and  the  maximum  pressure  variation  is  of  the  order  of  0.02 
atm.  In  the  far-held,  5L,  as  shown  in  Fig.  5b,  the  maximum  pressure  variation  is  of  the 
order  of  10“'*  atm.  and  is  reduced  by  a  factor  of  100  from  that  in  the  near-held.  In  Fig. 
5b,  the  difference  between  the  “decoupled  modeh’  and  the  “fully  coupled  model”  can  be 
obseved.  But  this  difference  is  of  the  order  of  10“^  times  the  pressure  variation  in  the 
near-held  and  is  within  the  accuracy  of  the  numerical  solution. 

From  a  computational  view  point,  it  is  important  to  compare  the  performance  of  the 
two  models  based  on  the  CPU  time  required  by  each  calculation.  In  the  linear  vibration 
regime,  the  “fvilly  coupled  model”  used  36000  seconds  of  CPU  time  on  a  Cray-ymp  to 
advance  the  calculation  by  10000  time-steps,  whereas  the  “decoupled  model”  used  only 
1000  seconds  for  the  same  calculation.  In  the  nonlinear  vibration  regime,  grid  rehnements 
were  needed  to  resolve  the  large  gradients  both  on  the  panel  and  in  the  radiation  held. 
Therefore,  in  order  to  advance  the  calculation  by  10000  time-steps,  72000  seconds  were 
used  by  the  “fully  coupled  model”  while  the  “decoupled  model”  used  only  1200  seconds. 

5.  Conclusions 

An  efficient  model  for  coupling  the  vibration  of  a  panel  to  the  on  surface  acoustic  radi¬ 
ation  was  derived.  The  model  uncouples  the  panel  vibrations  from  the  acoustic  wave  prop¬ 
agation  problem.  The  results  showed  that  this  model,  referred  to  a-s  “decoupled  model”, 
accurately  predicts  the  panel  response  and  acoustic  radiation  in  the  linear  and  nonlinear 
vibration  regimes  so  long  as  the  pressure  variation  in  the  how-held  remains  much  smaller 
than  the  ambient  pressure.  For  the  cases  studied  in  this  paper,  the  computational  cost  of 
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the  iiunieriral  integration  of  this  ino(h4  is  36  times  I'heaper  in  tlu'  linear  regime  and  GO 
times  cheaper  in  the  nonlinear  regime  than  the  cost  of  the  “fully  coujjh-d  model". 
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Figure  2:  Time  history  of  the  <lispla(  ement  of  the  i>auel  enter  for  an  excitation 
amplitude  of  100  dB  (linear). 
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Figure  5;  Time  history  of  the  racliated  i»ressure  along;  a  line  i)<'rpen<li<'nlai  to 
the  panel  renter,  (a)  1  inch  away  from  the  panel,  (h)  60  inches  away  from  the 
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